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We implement resonant single qubit operations on a semiconductor hybrid qubit hosted in a 
three-electron Si/SiGe double quantum dot structure. By resonantly modulating the double dot 
energy detuning and employing electron tunneling-based readout, we achieve fast (>100 MHz) Rabi 
oscillations and purely electrical manipulations of the three-electron spin states. We demonstrate 
universal single qubit gates using a Ramsey pulse sequence as well as microwave phase control, the 
latter of which shows control of an arbitrary rotation axis on the X-Y plane of the Bloch sphere. 
Quantum process tomography yields 7r rotation gate fidelities higher than 93 (96) percent around 
the X(Z) axis of the Bloch sphere. We further show that the implementation of dynamic decoupling 
sequences on the hybrid qubit enables coherence times longer than 150 ns. 


Isolated spins in semiconductors provide a promis¬ 
ing platform to explore quantum mechanical coherence 
and develop engineered quantum systems M. Silicon 
has attracted great interest as a host material for de¬ 
veloping spin qubits because of its weak spin-orbit cou¬ 
pling and hyperfine interaction, and several architectures 
based on gate defined quantum dots have been proposed 
and demonstrated experimentally PH eng. Recently, a 
quantum dot hybrid qubit formed by three electrons in 
double quantum dot was proposed m mi- and non- 
adiabatic pulsed-gate operation was implemented exper¬ 
imentally m, demonstrating simple and fast electrical 
manipulations of spin states with a promising ratio of co¬ 
herence time to manipulation time. However, the overall 
gate fidelity of the pulse-gated hybrid qubit is limited by 
relatively fast dephasing due to charge noise during one 
of the two required gate operations. Here we perform 
the first microwave-driven gate operations of a quantum 
dot hybrid qubit, avoiding entirely the regime in which 
it is most sensitive to charge noise. Resonant detuning 
modulation along with phase control of the microwaves 
enables a i r rotation time of less than 5 ns (50 ps) around 
X(Z)-axis with high fidelities > 93 (96) %. We also imple¬ 
ment Hahn echo [EMU an d Carr-Purcell (CP) [22] dy¬ 
namic decoupling sequences with which we demonstrate 
a coherence time of over 150 ns. We further discuss a 
pathway to improve gate fidelity to above 99%, exceed¬ 
ing the threshold for surface code based quantum error 
correction (23}. 

The quantum dot hybrid qubit combines desirable fea¬ 
tures of charge (fast manipulation) and spin (long co¬ 
herence time) qubits. The qubit states can be written 
as |0) = |^)| 5), where S denotes a singlet state in the 
right dot, and |1) = 1/a/ 3||)|T 0 ) — >/2/3|t)|X!_), where 
To and T_ are two of the triplet states in the right dot. 
The states |0) and |l) have the nearly same dependence 
on e in the range of detuning at which the qubit is op¬ 
erated (see Figs. lc,d), enabling quantum control that 
is largely insensitive to charge fluctuations. Moreover, 
electric fields couple to the qubit states and enable high 


speed manipulation on on nmm • Previously, we ex¬ 
perimentally demonstrated non-adiabatic quantum con¬ 
trol (DC-pulsed gating) of the hybrid qubit, where the 
manipulation and measurement scheme required the use 
of a detuning regime that is sensitive to charge noise 
(with £ near but not equal to zero—see Fig. Id) [H5] . 
Moreover, DC-gating requires abrupt changes in detun¬ 
ing. With a given bandwidth in the transmission line, 
pulse imperfections arising, e.g., from frequency depen¬ 
dent attenuation, lead to inaccurate control of rotation 
axes. In this work, we demonstrate resonant microwave- 
driven control and state-dependent tunneling readout of 
the qubit, which together overcome this limitation of DC- 
pulsed gating and enable full manipulation on the Bloch 
sphere at a single operating point in detuning that is 
well-protected from charge noise. 

The experiments here are performed in a double dot 
with a gate design as shown in Fig. la and with electron 
occupations as shown on the stability diagram of Fig. 
lb. The electron occupations and energy level alignments 
used for qubit initialization, readout, and microwave 
spectroscopy of the qubit states are shown schematically 
in Fig. lc. All the experiments reported here start with 
an initial dot occupation of (1,2) and with system in state 
|0), prepared at a detuning e ~ 230/ieV; this detuning is 
also used for measurement and corresponds to point M in 
Fig. lb. After initialization, we apply a microwave burst 
pattern at point O, which either coincides with point M 
or is reached through an adiabatic ramp in detuning (the 
latter case is illustrated in Fig. lb). Because of the pres¬ 
ence of the third electron in the pair of dots, the hybrid 
qubit states couples to time varying electric fields US], 
and qubit excitation or rotation occurs when the applied 
microwave frequency is resonant with the qubit energy 
level difference. The measurement point M is chosen so 
that the Fermi level of the right resorvior is in between 
the energies of |0) and 11), and we use the qubit state- 
dependent tunneling to project states |0) and |1) to the 
(1,2) and (1,1) charge states, respectively. Waiting at 
point M for « 10/is also resets the qubit to state |0), by 
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Figure 1. Microwave-driven coherent manipulation and readout of a hybrid qubit in a Si/SiGe double quantum 
dot device, a, SEM image and schematic labeling of a device lithographically identical to the one used in the experiment, b, 
Charge stability diagram near the (1,1)-(2,1)-(1,2) charge transition, showing the gate voltages used for microwave manipulation 
(O) and measurement (M). For clarity, linear background slope was removed from the raw charge sensing data, c, Schematic 
description of the qubit initialization, manipulation, readout, and reset processes, d, Energy E as a function of detuning e for 
the qubit states, calculated with Hamiltonian parameters measured in Ref. [18]. e, (inset) Probability P\ of the state to be |1) 
at the end of the driving sequence shown as a function of £ and the excitation frequency / of the microwave applied to gate R. 
In the main panel, the dashed green curve is the energy difference between the ground state and the lowest-energy excited state, 
as determined in mi f-j, Coherent Rabi oscillation measurements, f, P\ as a function of the voltage Vl and the microwave 
pulse duration tb with f— 11.52 GHz and excitation amplitude V ac = 400 mV. g, Linecut of Pi near Vl =-392 mV, showing « 
110 MHz Rabi oscillations. The red solid curve is a fit to an exponentially damped sine wave with best fit parameter T 2 * = 33 
ns. h-i, Rabi oscillation data with microwave amplitude 300mV (h) and 200mV (i). j, Rabi oscillation frequency /R a bi as a 
function of V ac with fixed / =11.52 GHz. The good agreement of a linear fit (red line) to the data is strong evidence that the 
measured oscillations are indeed Rabi oscillations, with the Rabi frequency proportional to the driving amplitude. 


tunneling an electron from the reservoir, if needed. Thus, 
the qubit state population following the microwave burst 
is measured by monitoring the current /qpc through the 
charge-sensing quantum point contact (Fig. la). Details 
of the measurement procedure and probability normal¬ 
ization are in Supplementary Information SI. 

We perform microwave spectroscopy of the qubit in¬ 
trinsic frequency—the energy difference 5E in Fig. Id— 
by applying the voltage pulse shown in the inset to 
Fig. le. The color plot in that figure shows the re¬ 
sulting probability of measuring state |1) after applying 
this pulse to initial state |0). The measured resonance 
and qubit energy dispersion agrees well with the green 
dashed curve, which is the calculated energy level dia¬ 
gram with Hamiltonian parameters measured in our pre¬ 
vious study [18]. As is clear from the color plot in Fig. le, 
the linewidth of the resonant peak narrows significantly 
at e > 200/ieV, becoming much narrower than the res¬ 
onance in the charge qubit regime (e 0) [28]. This 
linewidth narrowing corresponds to an increase in the 
inhomogeneous dephasing time, and it is this range in 
detuning that corresponds to the hybrid qubit regime. 


Applying microwave bursts to gate R in the hybrid 
qubit regime yields Rabi oscillations, as shown in Figs, lf- 
i. The Rabi frequency increases as a function of in¬ 
creasing microwave amplitude V ac (measured at the arbi¬ 
trary waveform generator), resulting in Rabi frequencies 
as high as 100 MHz. Fig. lj shows the power depen¬ 
dence of the qubit oscillations, revealing an oscillation 
frequency that is linear in the applied amplitude, as ex¬ 
pected for Rabi oscillations. The speed of the X-axis 
rotation demonstrated here is comparable to electrically 
manipulated spin rotations in InSb and InAs, which rely 
on strong spin-orbit coupling of the host material 1221130 ; 
here we achieve fast rotations solely through electric field 
coupling to the qubit states. This coupling is also highly 
tunable, since it is determined by the ground and ex¬ 
cited state inter-dot tunnel couplings m- Below, we 
characterize gates with the qubit frequency chosen to be 
« 11.52 GHz. 

We characterize the inhomogeneous dephasing time 
by performing a Ramsey fringe experiment, which also 
demonstrates Z-axis rotations on the qubit Bloch sphere. 
The microwave pulse sequence is shown schematically in 
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Figure 2. Ramsey fringes and two-axis control of the qubit, a, Schematic diagram of the pulse sequences used to 
perform universal control of the qubit. Both the delay t e and the phase 0 of the second microwave pulse are varied in the 
experiment, b-c, Experimental measurement of Z-axis rotation. In a Ramsey fringe (Z-axis rotation) measurement, the first 
X n /2 gate rotates the Bloch vector onto the X-Y plane, and the second X n / 2 gate (0=0) is delayed with respect to the first 
gate by t e , during which time the state evolves freely around the Z-axis of the Bloch sphere, b, Pi as a function of Vl and t e 
for states initialized near |Y). c, Pi as a function of t e with fixed Vl = —391.7 mV, showing « 11.5 GHz Ramsey fringes. The 
red solid curve is a damped sine wave with best fit parameter T 2 * = 11 ns. d, Effect of the phase 0 of the second microwave 
pulse on the state | Y) (by applying X^/ 2 on |0), black), and | —Y) (by applying X n / 2 and Z n on |0), red). The clear oscillation 
of Pi as a function of 0 in both cases demonstrates control over the second rotation axis by control of the phase 0. 


Fig. 2a. We first prepare the state | Y) = y/l/2(|0) + i|l)) 
by performing an X^/ 2 rotation. Z-axis rotation results 
from the evolution of a relative phase between states |0) 
and |1), given by ip = — t e SE/h , where t e is the time spent 
between the state preparation and measurement X ^/2 
pulses, the latter of which is used to project the Y-axis 
component onto the Z-axis. The final probability Pi is 
measured as described above. Fig. 2b shows the resulting 
quantum oscillations as a function of Vl, which controls 
the detuning energy, and t e . Fig. 2c shows a linecut taken 
near the optimal resonant condition (Vl ~ —392 mV), 
showing clear oscillations in Pi consistent with the qubit 
frequency of ~ 11.5 GHz. By fitting the oscillations to 
an exponentially damped sine wave (red solid curve), we 
extract an inhomogeneous dephasing time T 2 * = 11 ns, 
consistent with the value measured previously with non- 
adiabatic pulsed gating on the same device with similar 
intrinsic qubit frequency HU We estimate the typical 
tunneling-out time T 0 « 200 ns (see Supplementary In¬ 
formation SI), so that the inhomogeneous coherence time 
is not likely limited by electron tunneling to the reservoir 
during the measurement phase. 

Resonant microwave drive also enables arbitrary two- 
axis control on the X-Y plane of the Bloch sphere by 
varying the relative phase 0 of the X ^/2 pulses. Fig. 2d 
shows a measurement of Pi, demonstrating both two- 
axis control and phase control. Starting from a maxi¬ 
mum (minimum) Pi at 0 = 0, when we apply X 7r / 2 ^ 7 r /2 
(X n /2Z^Q^/2) on the state 10), Pi oscillates as a func¬ 


tion of the relative phase 0 that determines the axis of 
the ^tt/ 2,0 rotation. The deviation from an ideal sinu¬ 
soidal oscillation stems from limited phase resolution of 
our method of waveform generation (see Supplementary 
Information SI). 

We now turn to echo and dynamic decoupling pulse 
sequences. Fig. 3a shows Hahn echo mm and Carr- 
Purcell dynamic decoupling (CP) [22] pulse sequences. 
Provided that the source of dephasing fluctuates slowly 
on the timescale of the electron spin dynamics, insert¬ 
ing an X n pulse between state initialization and mea¬ 
surement, which is performed with X ^/2 gates, corrects 
for noise that arises during the time evolution. Fig. 3b 
shows a typical echo measurement. While keeping the 
total free evolution time r fixed at 10 ns, we sweep the 
position of the decoupling X n pulse to reveal an echo en¬ 
velope EU 122 ]. In Fig. 3b, the oscillations of Pi as a func¬ 
tion of St are at twice the Ramsey frequency (2/R amsey ~ 
23 GHz) and are well-fit by a gaussian decay (red solid 
curve). With r = 10 ns, we observe a maximum echo am¬ 
plitude ~ 0.3 near St ~ 0. The gaussian decay indicates 
that a significant component of the dephasing arises from 
low-frequency noise processes. 

Improvement in coherence times can be obtained by 
implementing CP sequences (see Fig. 3a), which use mul¬ 
tiple X n pulses inserted during the free evolution. Since 
the time scale for the CP sequence is typically longer 
than T 0 « 200 ns, an adiabatic detuning offset of am¬ 
plitude ~ 60 /ieV was applied during free evolution in 
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Figure 3. Increasing coherence time with dynamic decoupling sequences, a, Schematic pulse sequence for the 
measurement of Hahn spin echo and Carr-Purcell (CP) dynamic decoupling sequences that correct for noise that is static on 
the time scale of the pulse sequence DM2]- Note that for CP dynamic decoupling, free evolution is performed at ~ 60 fieV 
more positive detuning than the readout point, in order to prevent tunneling of the state |1) to the reservoir during the 
manipulation pulses and free evolution, b, Typical Hahn echo measurement with fixed total evolution time r — 10 ns, showing 
Pi as a function of Vl and delay time St of the pulse. The bottom panel shows line cut of Pi near Vl = -391.7 mV. The 
red solid curve shows a fit to a gaussian decay with fixed inhomogeneous coherence time T 2 * ~11 ns, obtained from the Ramsey 
fringe measurement of Fig. 2. The oscillations of Pi as a function of St at twice the Ramsey frequency (« 23.04 GHz) is clear 
evidence of spin echo, where we measure an echo amplitude of ~ 0.3 for r = 10 ns near St = 0. c, Pi as a function of r with 
fixed St = 0. The symbols shows the data, while the solid curves are the fit to the decay form Pi(r) — 0.5 + A e _ ^ T//T2 ^ with 
fixed exponent a = 2 for even numbers n of decoupling pulses ranging from 2 to 10. d, Coherence time T 2 as a function of 
n obtained from the fit of CP decay data to the decay form with fixed exponent a = 2 (black circles), 3 (orange triangles), and 
4 (purple diamonds). The resulting Th is insensitive to the choice of a to within the uncertainty. Applying the CP dynamic 
decoupling sequences increases the coherence time by more than an order of magnitude compared to T 2 « 11 ns. 


order to prevent electron tunneling to the reservoir. In 
the absence of dephasing, the CP sequence with an even 
number of X n pulses applied on the state |0) yields Pi = 
1. The measured Pi as a function of r, shown in Fig. 3c, 
decays exponentially due to dephasing: Pi(r) = 0.5 + A 
e -( r /? 2 )“, where a depends on the frequency spectrum 
of the dominant noise sources [33] . Fig. 3d shows the re¬ 
sults of fits as a function of the number n of decoupling 
X^ pulses with fixed a = 2, 3, and 4. The resulting co¬ 
herence time X 2 shows more than an order of magnitude 
improvement (> 150 ns) with n = 8, and the resulting 
times are approximately independent of the a used in 
the fit. Beyond n = 8 we typically observe a decrease in 
X 2 , which is likely due to accumulation of pulse imper¬ 
fections in the microwave or detuning pulses. We expect 
that optimization of microwave pulses can increase T 2 
further. 

We now present tomographic characterization of the 
microwave-driven hybrid qubit. To reconstruct the time 
evolution of the single qubit density matrix, we use the 
microwave pulse sequences shown schematically in Fig. 
4a to perform repeated state evolution under an X(Z) 


gate and perform independent X, Y, and Z axis projec¬ 
tive measurements. For state tomography under X(Z)- 
axis rotation, we prepare initial states near |0) and near 
| Y). After time evolution under the gate operation, we 
measure X, Y, and Z axis projections of the time-evolved 
Bloch vector using —Y(tt/2), X(tt/2), and identity oper¬ 
ations, respectively, and measure the resulting Pi. Note 
that the pulse sequences shown in Fig. 4a represent state 
tomography in the rotating (laboratory) frame for Rabi 
(Ramsey) oscillations, because the phase of the second 
7r/2 pulses for the Rabi oscillation tomography evolves 
as the length of Rabi manipulation % is increased, while 
the second microwave pulse in the state tomography of 
the Ramsey fringes has fixed relative phase with respect 
to the first microwave pulse. For comparison, we also 
present the state tomography of the Rabi oscillation per¬ 
formed in the lab frame in Supplementary Information 
S2. Figs. 4b and 4c show X (black circles), Y (green tri¬ 
angles), and Z (orange squares) axis projections of the 
time-evolved Bloch vector under continuous X (b) and Z 
(c) axis rotation gates. 

Based on the density matrices obtained from the state 
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Figure 4. State and quantum process tomography (QPT) of the ac-driven hybrid qubit, a, Schematic of microwave 
pulse sequences used for the tomographic characterization of continuous Bloch vector evolution under X (left panel) and Z (right 
panel) axis rotations, b-c, X-axis projection Px (green), Y-axis projection Py (orange), and Z-axis projection Pz (black) of 
the Bloch vector evolution under X (b) and Z (c) axis rotations performed in rotating (lab) frame (see main text for more 
discussion), d-g, Real and imaginary parts of the process matrices x EU i n the Pauli basis {I, A, Y, Z} for processes X n / 2 
(d), X n (e), Z n / 2 (f), and Z n (g) obtained by quantum process tomography with maximum-likelihood estimation 18 34l (35] 
(real: open circles, imaginary: open triangles) compared to corresponding ideal processes (bars). The error bars represent the 
standard deviation of the result obtained by 10 distinct input and output density matrices chosen from the state tomography 
data. Process fidelities defined by P p = Tr(xideaiXML-QPT 1 ) are 93% for X n and 96% for Z n processes, respectively (left panels 
of e and g). 


tomography, we implement quantum process tomography 
(QPT) to extract fidelities of single qubit gates on the ac- 
driven hybrid qubit through the relation PI 1531155], 

4 

£ (P) = X E m pElxmn, (1) 

m,n= 1 

where £(p) is the density matrix specifying the output for 
a given input density matrix p, the E m are the basis op¬ 
erators in the space of 2 x 2 matrices, and x is the process 
matrix. Experimentally, four linearly independent input 
and output states are chosen from continuous evolution 
of the state under X and Z axis rotation available from 
the state tomography data set, and the maximum likeli¬ 
hood method nails] is used to determine Figs. 4d-g 
show the results of QPT (symbols) performed on the tt/2 
and 7r rotations around the X and Z axes and comparison 
to corresponding ideal rotation process matrices (bars). 
The error bars of length ~ 0.01 to 0.02 represent the 
standard deviation of the experimental result obtained 
by 10 distinct input and output density matrices chosen 
from the state tomography data. The process matrices 
X obtained from QPT in the Pauli basis {I, cr x , cr y , cr z } 
yield process fidelities F p = Tr (xideaix) of 93% and 96% 


for 7r rotations around X and Z axes, respectively. Com¬ 
paring these results to the process fidelities of 85% and 
94% for X and Z-axis rotation reported previously for 
the non-adiabatic DC-pulse gated hybrid qubit [18], we 
find more than a factor of two reduction in the X-axis 
rotation infidelity. 

The improvement in overall fidelity of the AC-gated 
quantum dot hybrid qubit demonstrated here compared 
to DC-pulsed gating stems mainly from (1) elimination 
of the need to enter the regime in which the qubit is 
sensitive to charge noise by using resonant manipulation 
and tunneling-based readout, and (2) reduced rotation 
axis and angle errors because resonant driving with fixed 
frequency enables more accurate control of these quan¬ 
tities. The AC driving in this work was performed by 
resonantly modulating the energy detuning between the 
dots. For this type of modulation, the ratio of manip¬ 
ulation time (Rabi period) to coherence time depends 
strongly on the strength of ground and excited state 
tunnel couplings m- Thus we expect that further fi¬ 
delity improvement can be achieved by adjusting tunnel 
couplings. Moreover, recent theoretical work suggests 
that dynamically modulating tunnel coupling instead of 
detuning can enable Rabi frequencies exceeding 1 GHz 
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while keeping long coherence times, enabling achievement 
of gate fidelites exceeding 99% [36 . Finally, implemen¬ 
tation of dynamic decoupling sequences can be useful to 
understand the noise spectrum of the system as well as 
to further enhance the coherence time. 

Methods 

The details of the Si/SiGe double quantum dot de¬ 
vice are presented in Refs. [32, 37. We work in the 
region of the charge stability diagram where the va¬ 
lence electron occupation of the double dot is (1,1) or 
(1,2), as confirmed by magnetospectroscopy measure¬ 
ments mami- All manipulation sequences, including the 
microwave bursts, are generated by a Tektronix 70002A 
arbitrary waveform generator and are added to the dot¬ 
defining dc voltage through a bias tee (Picosecond Pulse- 
labs 5546-107) before being applied to gate R. We map 
the state |0) and |1) to (1,2) and (1,1) charge occupa¬ 
tion states respectively, leading to conductance changes 
through the quantum point contact (QPC). We measure 
with a lock-in amplifier (EG&G model 7265) the differ¬ 
ence in conductance with and without the applied mi¬ 
crowave burst. When converting time averaged conduc¬ 
tance differences to the reported probabilities, tunneling 
between the (1,2) and (1,1) charge states during the mea¬ 
surement phase is taken into account using the measured 
times for tunneling out of (T 0 ~ 200 ns) and into ( T\ ~ 2.1 
/is) the dot. Supplementary Information SI presents the 
details of the measurement technique and the probability 
normalization. 
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Supplementary Fig. SI. Measurement of state popu¬ 
lations using energy dependent loading and unload¬ 
ing. a, Schematic charge stability diagram as a function 
of gate voltages Vl and Vrp. The close-up view near the 
(1,2) to (1,1) charge transition shows a region (orange box) 
where state manipulation using resonant microwaves as well 
as robust readout using energy dependent loading and un¬ 
loading can be performed. In this readout window, applied 
microwaves are resonant with the qubit energy spacing; in ad¬ 
dition, the Fermi level of the right reservoir lies between the 
qubit states, so that qubit state |1) leads (after tunneling) to 
a QPC current level corresponding to the (1,1) charge state, 
whereas state |0) leads to that of (1,2). b, Experimentally 
measured variation of the readout window as a function of 
excitation microwave frequency / with a fixed pulse duration 
of « 5 ns. The length in gate voltage space of the readout 
window increases at higher excitation frequency, reflecting the 
saturation of the qubit level spacing at large detuning, c, 
Schematic pulse sequence to manipulate the qubit states and 
measure the population of state |1). Note that the adiabatic 
in- and out- ramps do not induce state population change; 
this baseline change is included even without microwaves, in 
order to reduce the background lock-in signal that would oth¬ 
erwise arise from capacitive crosstalk, d, Measurement of 
uncalibrated probability of state 11) , Pi, m , as a function of 
the measurement time £ m . The solid red curve shows a fit 
to the solution of rate equations Eqs. S1-S7 with excitation 
efficiency a = 0.93, tunnel-out time T 0 « 200 ns, and reset 
time Ti « 2.1 jus. 


SUPPLEMENTARY INFORMATION 

SI. DETAILS OF READOUT AND 
CONVERSION TO STATE PROBABILITIES 

We use the qubit state-dependent loading and unload¬ 
ing rates to implement robust readout of the state |1) 
probability. Supplementary Fig. [ST} i shows a schematic 
charge stability diagram (left panel) and a magnified re¬ 
gion near the (1,1)-(1,2) charge transition (right panel). 


The green dashed region (tunneling condition) denotes 
where the Fermi level of the right reservoir lies between 
the qubit energy levels. In this region, the energy of state 
|1) is above the Fermi level of the resorvoir, so in this 
state an electron tunnels out to the right reservoir with 
characteristic tunneling time T 0 , leading to an Iqpc level 
corresponding to the (1,1) charge configuration, whereas 
state |0) remains in the (1,2) charge state. In our exper¬ 
iment, the qubit energy splitting is typically >11 GHz, 













7 


at least four times larger than the thermal broadening of 
the reservoir Fermi level (^140 mK or « 2.8 GHz), which 
enables high signal contrast between the qubit states. 
Moreover, it has been shown that the qubit relaxation 
time Ti from |1) to |0) is > 100ms [16]; thus, the read¬ 
out fidelity and overall coherence time is not limited by 
Ti. Excitation from |0) to |1) is performed by applying 
fixed-frequency microwaves to gate R. The excitation is 
possible when the detuning is such that the microwaves 
are resonant with the qubit frequency, which we denote 
as a purple strip (resonant condition) in Fig. Sir. The 
intersection of these conditions (the read out window, see 
the orange box in Fig. SI) defines the region of the stabil¬ 
ity diagram where the applied microwaves are resonant 
with the qubit level splitting and also the final state can 
be read out. Supplementary Fig. |S1 |d shows the variation 
of the read out window identified by sweeping the volt¬ 
ages on gates L and RP while applying microwave bursts 
on gate R. 


Manipulation sequences, including the microwave 
bursts, are generated using a Tektronix AWG70002A 
arbitrary waveform generator (AWG) with a maximum 
sample rate of 25 Gs/s and an analog bandwidth of about 
13 GHz. To generate high frequency (> 10 GHz) mi¬ 
crowave signals, we adopted two approaches: (1) For 
experiments that do not require relative phase control 
between sequential microwave pulses, we adjust the sam¬ 
ple rate of the AWG to generate a sequence of triangular 
wave segments so that the resulting sequence has a fun¬ 
damental frequency matching the desired excitation fre¬ 
quency. Higher harmonic components are filtered using 
a low pass filter with a cutoff frequency of about 18 GHz. 
Additionally, transmission lines from room temperature 
electronics to the sample also provide effective attenua¬ 
tion for filtering out > 20 GHz harmonic components. (2) 
For relative phase control, e.g., the data in Fig. 2d, we 
increased the effective sampling rate by combining two 
channels of the AWG, with the second channel output 
delayed by one half of the sampling time. By construct¬ 
ing odd (even) numbers of sampling points using the first 
(second) channel of the AWG, we can directly synthesize 
high frequency >10 GHz waveform limited only by the 
analog band width (13 GHz) of the instrument. 

For measuring changes in qubit state occupation prob¬ 
abilities resulting from fast microwave bursts, we use 
the general approach described schematically in Fig.|Slfc, 
where we measure the difference between the QPC con¬ 
ductance with and without the manipulation pulses. The 
data are acquired using a lock-in amplifier with a refer¬ 
ence signal corresponding to the presence and absence of 
the pulses (lock-in frequency « 777 Hz). We compare the 
measured signal level with the corresponding (1,1)-(1, 2) 
charge transition signal level, calibrated by sweeping gate 
L and applying square square pulses with frequency of 
777 Hz to gate L. 


To calibrate the state probabilities, we account for the 
two time scales that affect the overall voltage signal dur¬ 
ing the measurement phase: (1) the tunneling time T 0 


from |1) to (1,1), which determines the onset of the /qpc 
change, and (2) the tunneling time T\ from (1,1) to jo), 
which determines the effective reset time. As shown 
below, the previously measured [16 Ti-time of ~ 140 
ms is much longer then T 0 and Tp thus, the effect of 
T\ during the measurement phase is negligible. To ex¬ 
tract these time scales from experimental measurements 
at the readout window, as discussed below, we consider 
the three-state rate equation describing the occupations 
no(£), ni(t), and n( 1?1 )(£) of the corresponding states |0), 
|1), and (1,1): 

d n 0 {t)_ 

— ^(1,1 )v')/Ti , (SI) 

rfn(1 ^ )(t) = — n (i,i)(t)/Ti + n\{t)/T 0 , (S2) 

^ = (S3) 

The boundary conditions appropriate to repeated pulses 

with duration t m are 

M°) + m(0) + n (1}1 )(0) = 1 , (S4) 

^(1,1) (0) = n (M )(£m) , (S5) 

ni(0) = cm 0 (£ m ) + (1 - a)ni(t m ) , (S6) 

n 0 (0) = cmi(£ m ) + (1 - a)n 0 (t m ) . (S7) 

The rate equations describe relaxation of 11) to (1,1) with 
tunneling time T 0 and relaxation of (1,1) to |0) with tun¬ 
neling time Ti ; the boundary conditions follow because 
the microwave pulses convert |0) to |1) with X n gate fi¬ 
delity a and vice versa, and do not affect the occupation 
of the (1,1) configuration. 

Experimentally, we measure the uncalibrated proba¬ 
bility Pq m of state |1) after an X(7r) pulse, which will 
depend on the measurement time t m (see Supplemen¬ 
tary Fig. |Sl]y), and compare the data to the solution of 
the three-state rate equations. Supplementary Fig. [STJi 
shows Pi : m as a function of £ m , showing the expected re¬ 
duction of signal size arising from the reset to state |0). 
Since the value of a {X n gate fidelity) is not known a 
priori , we fit the data to the solution of the rate Eqns. 
S1-S7 with initial guess of <a=0.95, and perform the pro¬ 
cess tomography (QPT) to extract X n gate fidelity. We 
repeat the fit procedure to self-consistently determine a 
from the normalization fit and X n gate fidelity from the 
QPT, with a tolerance in the difference of =b 0.01, and 
obtain a converged value of <a=0.93 after two iterations, 
consistent with the X n gate fidelity in the main text. 
The red solid line of Supplementary Fig. |Sl] i shows the 
converged fit, which yields T 0 « 200 ns and T x 2.1 /i s. 
The probability data in the main text use t m ~ 11 fi s, 
long enough to ensure state reset after the manipulation 
sequences. Using these time scales allows us to correct 
for state reset by appropriate normalization of the data. 
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Supplementary Fig. S2. State tomography of the Rabi 
oscillation in the laboratory frame, a, Schematic Bloch 
sphere representation of the expected state evolution upon 
resonant microwave drive in the laboratory frame, b, Mi¬ 
crowave pulse sequences used to measure the X-axis projec¬ 
tion of the Bloch vector during the Rabi oscillations. Note 
that the measurement pulse {—Y^/ 2 ) has a fixed phase shift 
of -90 degrees with respect to the start of the Rabi drive pulse, 
which has length tb- c, X-axis (red circles) and Z-axis (black 
triangles) projection of the Bloch vector during a Rabi pulse 
obtained by lab frame tomography measurements as a func¬ 
tion of tb- The data is shown as a function of tb up to a 
time corresponding to an X n rotation. The X-axis projec¬ 
tion is clearly consistent with a spiral trajectory of the Rabi 
nutation in the lab frame, which is shown schematically in a. 


S2. STATE TOMOGRAPHY OF RABI 
OSCILLATIONS IN THE LABORATORY FRAME 


The state tomography of Rabi oscillations in the main 
text are performed in the rotating frame. Here we also 
show tomographic measurements in the lab frame. Sup¬ 
plementary Fig.|S2^t shows a schematic of the evolution of 
the Bloch vector upon resonant microwave excitation. To 
implement measurement in the lab frame after microwave 
manipulation, we fix the phase of the measurement ir/2 
pulse and apply it immediately after a Rabi drive pulse 
of length tb- Supplementary Fig. |S2 | d shows, as an ex¬ 
ample, the microwave sequence for the X-axis projection 
measurement of the continuous Rabi oscillation. Note 
that Larmor phase accumulation does not affect the Z- 
axis projection measurement, so that the rotating and 
lab frame measurements of the Z-axis component of the 
Bloch vector are identical. The result of the lab frame 


state tomography is shown in Supplementary Fig. S2 


The rapid oscillation of the X-axis projection P x of the 
probability varies at the qubit frequency of « 11.52 GHz, 
starts and ends at probability near 1/2, and is consis¬ 
tent with the expected time evolution of the qubit state 
under resonant microwave drive, showing good control 
over the measurement axis in this experiment. More¬ 
over, measurement of state evolution in the lab frame can 
be a useful tool to understand decoherence of a driven 
qubit ma so]. The resolution (i.e., the point density) 
of the measurement of this rapid precession is currently 
limited by the sampling time of the arbitrary waveform 
generator (Tektronix 70002A, set to 43.4 ps). 
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